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Q^' ABSTRACT 
<' 

I MeV seed photons produced in shocks in a variable ultra-relativistic outflow gain 

energy by the Fermi mechanism, because the photons Compton scatter off relativistically 
^ colliding shells. The Fermi- modified high-energy photon spectrum has a non- universal 

\^ • slope and a universal cutoff. A significant increase in the total radiative efficiency is 

. possible. In some gamma ray bursts, most of the power might be emitted at the high- 

. energy cutoff for this mechanism, which would be close to 100 MeV for outflows with a 

, mean bulk Lorentz factor of 100. 

O 
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, 1. Introduction 



The most common model of gamma-ray burst (GRB) sources involves a relativistic outflow 
in which shocks occur and radiate away a fraction of the bulk kinetic energy. For typical model 
parameters the synchrotron emission peaks around 1 keV in the comoving frame, or ~ 100 keV in 
the observer frame. Generally, this is considered to be the primary spectrum observed. However, 
in internal shocks in the neighborhood of the flow photosphere, and also in external shocks in 
some cases (e.g. Madau, Blandford, & Rees 2000), the shocks can have a non-negligible Thomson 
scattering depth, which results in upscattering of these primary photons. Single scattering on 
individual shock-accelerated electrons with Lorentz factors 7e ~ 300 produces photons with energies 
~ 7g keV in the comoving, or Tj'^ ~ 10r2 GeV in the observer frame, where F = IOOF2 is the 
bulk Lorentz factor. 

Here we concentrate on a different, multiple scattering component. This is associated with 
mildly relativistic motions of different ejecta shells or turbulent cells resulting from multiple shock 
interactions, which contain more energy than the shock-accelerated highly relativistic electrons. 
Multiple interacting shells are naturally expected in internal shocks, and also in external shocks 
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when a longer lasting modulated outflow runs into the flrst decelerated shell (e.g. Fenimore & 
Ramirez-Ruiz 2000; Kumar & Piran 2000). 

Repeated scatterings using the energy of these bulk motions boost the photon energy through 
the equivalent of the Fermi acceleration mechanism of particles (Blandford & Payne 1981). This is 
related to Thompson's (1994) photon scattering off Alvfen waves, but our mechanism relies instead 
on relative bulk motions. It differs also in using synchrotron photons instead of thermal photons 
as its source term, and hence leads to different characteristic energies. This bulk Comptonization 
results in a spectrum extending at least up to ^ MeV in the comoving frame and ^ 100r2 MeV in 
the observer frame. The spectral power or luminosity per decade can increase as steeply as linearly 
in the photon energy. This provides a natural explanation for those GRB spectra (e.g. Preece et 
al 1999) which show a positive i^Fu slope above the MeV range (generally < +1), which cannot be 
explained by direct synchrotron radiation from Fermi shock-accelerated electrons. 

The component made up of bulk-scattcrcd photons can extend up to a maximum observed 
energy ~ 100r2 MeV. Beyond this energy, Klein-Nishina and electron recoil effects set in, and the 
spectrum reverts to being dominated by the seed spectrum (with negative or flat power law slope) 
of the unscattered photons above the synchrotron peak. 

In this work we adopt a test photon model, which assumes no back reaction on the plasma. 
Several potentially important effects are neglected (upscattered photons can heat electrons in the 
colliding shells and produce pairs, light pressure ensures that the total comoving energy of the 
scattered photons cannot exceed the total kinetic energy of the relative shell motions). These 
effects will be investigated elswehere (Gruzinov, Mcszaros , &: Rees 2000). The simplified approach 
that we use here retains the essential properties of the bulk motion comptonization phenomenon, 
and allows us to explore the main qualitative features it introduces in the spectra. 

In §2 we specify the GRB model, in §3 we give an analytical model of Fermi acceleration of 
photons by colliding shell, and in §4 we describe our Monte Carlo simulations. The results are 
discussed and related to current and future observations in S5. 



2. The GRB internal shock model 

As a specific example to illustrate the effect, we restrict ourselves here to the internal shock 
model. We use a standard set of parameters for the ultra-relativistic outflow model of GRBs (e.g. 

Mcszaros & Rccs 2000), i.e. luminosity L = 10^^L52 erg/s, L52 ~ 1; terminal Lorcntz factor 
r ~ 100; variability time scale = W~'^tms s, with tms ~ 1- With these parameters the shock 
between shells ejected at typical time intervals ty with AF F occurs at a radius R ~ T^cty. To 
determine the optical depth to Thomson scattering, we estimate the proper density of the wind as 



47ri?^F^m„c3 
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The proper radial width of colhding shehs is ~ i?/r, and the optical depth is 

T ~ naxR/T ~ - — I -f^ I ■ (2) 

We see that optical depths of ^ 0.1 are natural. In fact, values of r ~ 1 are not purely coincidental, 
since in the standard GRB model the parameters are chosen so as to make r < 1 enabling a non- 
thermal spectrum, and the "preferred" model uses smaller values of T (which do not require overly 
low baryon loads), thus making r ^ 1. The observational fits to the models provide reasonable 
support for such a choice of parameters (e.g. van Paradijs, Kouveliotou & Wijers, 2000). 

Most of the seed photons are emitted at the synchrotron peak. Assuming mildly relativistic 
internal shocks, and magnetic field and electron energies equal to a fraction and of their 
equipartition values, the synchrotron peak frequency in the comoving frame is 



3. Analytical model of photon acceleration by colliding shells 



The basic features of the bulk Comptonization process can be understood by means of a simple 
analytical model, which reproduces the essence of the Monte Carlo simulation discussed in §4. This 
model is sufficient to show that the high-energy slope of i^F^ can be positive, and gives a value for 
the cut off energy. 



observer frame 



comoving frame 







Fig. 1. — A photon scatters off relativistically colliding shells. As seen by shell A, the scattered 
photon has a small angle of incidence. 

Consider a photon between two shells of equal optical depth r <C 1 (Figure Q). The shells 
collide with an ultra-relativistic relative Lorentz factor 7 3> 1 Q. Initially, in the frame comoving 
with shell B, the photon frequency is v and the incidence angle is 0. We assume the Thomson 
regime, hu <C mgC^. The probability that the photon passes through shell B without scattering is 



^7 = i (^7^ + Y^)- The case of 7 ~ 1 requires a more cumbersome analysis. At our accuracy, we just assume 
that the ultra-relativistic solution is a useful approximation in the mildly relativistic regime. 
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1 — r. The photon scatters forward with probabihty t/2, and it scatters backward with probabiHty 
t/2. In the frame comoving with shell A, the backward scattered photon has a small incidence 
angle. The average (over the scattering kernel) of the Lorentz-transformed frequency is 

= j:^,d9sin9{l + cos-e) 

Thus, after one scattering the total energy of photons is changed by a factor of (25/32)r7, and 
the frequency is changed by a factor of (25/16)7. Then the observed spectrum is a power law with 
the luminosity per frequency octave uF,^ oc v^, with 

^ log(25T7/32) 

log(257/16) ■ ^ ' 

The slope is positive if r7 > 1.28. The steepest slope is /? = 1, which is achieved in the limit 
7 ^ T~^. This model is similar to that for the ultrarelativistic isotropic multiple scattering (c.f. 
Rybicki & Lightman 1979, p. 212) where the slope is /3 = 1 + log(T')/log(7'^). In the last equation 
we should use r' ~ r and 7' = ^^(1 + 7) /2, which is the Lorentz factor of the colliding shells in the 
center of mass frame. 

A finite convergence time of the two geometrically thin shells does not introduce a limit on the 
maximum photon energy, as can be seen from the equivalent of Zeno's paradox (since the photons 
always travel faster than the shells approach, the number of photon reflections can in principle grow 
arbitrarily large for a shrinking separation). However, the effective number of scatterings should be 
limited, introducing a corresponding spectral cut off at high energies, due to (i) a gradual Klein- 
Nishina decline of the cross section; (ii) a sharp cut-off of the energy of scattered photons due to 
electron recoil (the maximal energy of a photon backward scattered from shell B observed in frame 
A is 7meC^); (iii) radiation pressure back reaction effects (if /3 > 0, the total energy of the scattered 
photons builds up and the photon pressure can prevent the collision of the shells). 



4. Monte Carlo simulations 

Monte Carlo simulations of the photon acceleration model sketched in Figure ^ were performed 
in order to check the simple estimates of §3. These confirm that the analytical results, including 
equation (|5|), are a useful approximation in the mildly relativistic regime. They also allow us to 
find in more detail the shape of the Compton recoil / Klein-Nishina cut off (see Figure |^. The 
Monte Carlo simulations were done using two approximations, which simplify the problem without 
affecting significantly the result: polarization effects are neglected, and multiple scattering during 
one passage through a shell is neglected. 

Initially, a seed synchrotron spectrum is released. The seed spectrum is isotropic in the center 
of mass frame, rising at low energies (/3 = 4/3), with a smooth transition to a spectrum decreasing 
at high energies (where we take as an example /? = —0.5). 
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Fig. 2. — Monte Carlo simulations of bulk Comptonization between two converging slabs (relative 
Lorentz factor 7 = 5), showing the comoving flux per logarithmic photon energy interval for three 
different optical depths r. The observed spectrum is the same multiplied by the mean bulk Lorentz 
factor r = lO^Pa. 

In Figure |2| we show a spectral calculation for two shells similar to those of Figure using 
a relative Lorentz factor 7 = 5. The unscattered seed synchrotron spectrum is marked "r = 0". 
The spectrum of all photons that finally escape outwards (to the right) through the shell B, as 
measured in the (comoving) frame B, is shown in Figure |^ for different values of r. The Compton 
recoil / Klein-Nishina effects are noticeable at about 1 MeV comoving energies, and there is a sharp 
Compton recoil cut off at hv / {rrieC^) = 7. In the observer frame this spectrum would be blueshifted 
by the factor F^. 

The dependence on the relative Lorentz factor is illustrated in Figure |^, showing the effect for 
a fiducial relative Lorentz factor 7 = 2, which may also be characteristic of the reverse shock in an 
external deceleration shock scenario. As expected, the steepening remains present but becomes less 
strong as weaker shocks are considered. The cut-off, however, is inherent to the scattering physics, 
and although it has a tendency to decrease somewhat with decreasing 7, its order of magnitude 
remains in the neighborhood of hu ~ mgC^ ~ 1 MeV (comoving frame). 



-3 



-2 



-1 

log(hi//m^c^) 







Fig. 3. — Same as figure |2|, 



comoving frame spectrum for a case with 7 = 2. 
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Discussion 



Tlie Fermi acceleration of photons in the standard GRB fireball model has several consequences 
of theoretical and observational interest. The first is that it can naturally produce vFi, spectra 
which are harder than the input synchrotron spectrum, including the possibility of an increasing 
i^Fi, above the usual break found in the Band parameterization of spectra (Band et al 1999, Preece 
et al 1999). The latter is a property that is hard to obtain with a synchrotron model. Such rising 
uFiy oc are expected, in this model, to have /5 ^ 1, and this implies that in some GRB most of 
the energy is at energies hi' ~ 100(r/102) MeV, well above the BATSE instrument band on the 
Compton Gamma Ray Observatory (CGRO). The current observational situation is that BATSE 
finds approximately 16 % of the spectra to be rising at ^ 1 MeV, and the observed rise is not faster 
than f3 = 1 (Preece et al 1999). These fits generally cut off above 1.8 MeV, and the break is usually 
not much lower than this, so there is some uncertainty. The COMPTEL instrument on CGRO, 
sensitive up to 30 MeV, has analyzed ^ 30 bursts (Schoenfelder et al 2000), and an analysis of the 
slopes indicates in several cases z^Fj, slopes /? ~ 0, with one burst of /3 ~ 0.5 (Kippen et al 1999). 
The EGRET experiment on CGRO has detected ~ 30 bursts with the scintillation counters in the 
1-200 MeV range, and ~ 7 bursts with the spark chambers in the 100 MeV-30 GeV range. In 
this range the spectra are largely noise-dominated (Schaefer et al 1998, Bromm & Schaefer 1999). 
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However, the scintillation spectral slopes (Catelli, Dingus & Schneid, 1997) are compatible with 
P ^ 0, although some could be positive and others negative, which is compatible with the analysis 
of spark chamber data (e.g. Sommers, 1994; Hurley et al, 1994). Large area detectors such as 
GLAST should be able to obtain more definite answers in the 20 MeV-300 GeV range. 

Another implication of the photon acceleration described here is that it provides a natural 
mechanism to increase the efficiency of conversion of baryon bulk motion into photon energy. This 
is of interest since in general the internal shock synchrotron efficiency for radiating in the BATSE 
band is limited to 1-10% (Kumar 2000; Spada, Panaitescu & Meszaros 2000; see however Fenimore &: 
Ramirez- Ruiz 2000, and observation-based estimates by Freedman & Waxman 2000). The increase 
in the radiative efficiency is simply given by the increase in the value of vFi, at different energies, 
or by the integral / F^diy in the range of interest. In the generic examples shown, this increase is 
substantial. Of course, the spectra shown in Figures ^ and ^ are test photon spectra, which do 
not take into account the back-reaction of radiation. The latter can become important when a 
substantial fraction of the bulk energy has been converted to radiation through Fermi acceleration, 
and the natural limit for this effect can be estimated as ~ 50% radiative efficiency. 

The calculations presented here are meant to illustrate the consequences of photon acceleration 
by bulk motions. A self-consistent calculation is needed in order to explore the back-reaction of the 
radiation pressure and pair formation on the shell dynamics. Pair formation has an angle averaged 
cross section a^^ ~ (Tt/8 at a threshold which is similar to that for Klein-Nishina effects. It is not 
very important in low compactness situations (shocks at large radii), while in high compactness 
cases the exponential tail of the photon distribution would lead to a pair cascade, which lowers the 
cut off in the comoving spectra of Figures |,| to ^ 0.5 MeV, possibly with some pile-up of photons 
at this energy (Gruzinov, Meszaros & Rees, 2000). 
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NASA NAG5-2857, the Guggenheim Foundation and the Institute for Advanced Study. 
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